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Abstract 

The normal state optical conductivity is calculated for a layered metal with 

two layers per unit cell coupled through a transverse hopping matrix element 

t±. The optical response involves an interband term in addition to the more 

familiar intraband term which leads to the usual Drude form. The interband 

term is only weakly temperature dependent, even for an inelastic scattering 

rate which is linear in T. It gives a c-axis response which extends in frequency 

over the entire band width although there can be structure on this energy 

scale which reflects details of the electronic structure. In particular, at low 

energy, the c-axis response can develop a gap or pseudogap as the temperature 

is lowered. At high temperature, a Drude response will be seen only if the 

intraband transitions, which are of order tj_, become important compared 

with the interband transitions which are of order t±. 
74.25.Nf,74.25.Jb7,74.72.-h 
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I. INTRODUCTION 



Models of the high T c oxides often start with a single isolated Cu02 layer. Other struc- 
tural elements within the unit cell are usually labelled as charge reservoirs or barrier layers 
and are ignored. Yet it is precisely the nature of the coupling between the different layers 
that determines c-axis transport properties, which are found to display a rich variety of 
behaviours and are often anomalous.!!! The anomalous nature of the c-axis conductivity has 
been variously interpreted as suggesting that interlayer coupling is an essential piece of the 
superconducting mechanism!, that the interlayer coupling is incoherent due to impurity or 
phononic scattering!^ or thermal fluctuations!, or that the Cu02 layers are in a non-Fermi 
liquid stated. In this work, we suggest that the c-axis optical conductivity in YE^CuaOz-^ 
(YBC0 7 _,5) can be explained by proper consideration of the multilayered structure of the 
unit cell. 

In considering c-axis properties, it is important to distinguish between the coupling 
between various conducting layers within a unit cell (which can contain several planes) 
and the intercell coupling which could involve some barrier layer. It is this latter coupling 
that probably governs the size of the anisotropy observed between c- axi s0@'0HHI and afr-plane 
properties. On the other hand, the large in-plane anisotropy between a- and 6-directionsEHl3 
(along the CuO chains) observed in YBCO is related more closely to the properties within a 
unit cell. The actual situation is clearly quite complex. For example, the unit cell in YBCO 
consists of a bilayer of two Cu02 planes separated by a CuO chain layer. In addition, 
the chains are only completed in YBa2Cus07 and the effect of oxygen doping on the chain 
Fermi surface is not well understood nor is the partition of holes between planes and chains. 
Because of these uncertainties, it is necessary, at this stage, to use a simplified model and 
to set specific but limited aims. 

Having recognized that several distinct transverse hopping matrix elements come into a 
complete description of the c-axis properties of the oxides, we will, nevertheless, limit our- 
selves here to a model of two layers per unit cell coupled through a single transverse matrix 



2 



element t±. With YBCO in mind, one of the two layers will be assumed to have tetragonal 
symmetry and model a Cu02 plane, while the other will be taken to have orthorhombic sym- 
metry and represent a CuO chain. While this model is admittedly crude, it does allow us 
to examine the role of interband transitions on the optical conductivity. We are interested 
in addressing two questions: what is the magnitude of the interband contribution to the 
conductivity compared to the intraband (or Drude) contribtuion, and how different is the 
frequency dependence of the interband contribution from that of the Drude contribution? 

Not surprisingly, we find that the interband contributions are of order t\ and are therefore 
relatively unimportant for the a- and 6-axis optical conductivities. On the other hand, the 
results are reversed for the c-axis conductivity: the interband contributions are of order 
t\ and dominate the intraband contributions — which are of order t\_ — for weak interlayer 
couplings. It is not surprising, then, that the c-axis conductivity should have a non-Drude 
frequency dependence. In the work which follows, we will examine this frequency dependence 
and compare it with experiment. 

The paper is structured as follows. In section [TT], general expressions for the optical 
conductivity of a system with two layers per unit cell are derived. In section |TJ a specific 
model (the plane-chain model) which is suitable for YBCO is introduced. Numerical results 
are given for the conductivity which we present separately for a- b- (along the chains) and c- 
(perpendicular to the planes) directions. In section |TV], the expression for the conductivity 
derived in the previous section is reduced analytically, with the intention of highlighting 
the two types of contribution (interband and intraband) to the conductivity. One of the 
important results of this section is to show how the different contributions to the conductivity 
depend on the chain-plane coupling t±. In section [V|, expressions for the conductivity in the 
case of a bilayer — consisting of two identical but unevenly space planes in each unit cell — 
are derived. The calculation is interesting because YBCO, as well as many other cuprate 
superconductors, contains a Cu02 bilayer in the unit cell. One of our main conclusions in 
this section is that the coupling between the Cu02 planes is not as likely to be the source 
of the broad background seen in the c-axis optical conductivity as is the coupling between 
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the C11O2 planes and the chains. Section [V| consists of a short discussion of sum rules. A 
long conclusion, which includes some further discussion and a summary, is to be found in 
the final section. 



II. CONDUCTIVITY IN A LAYERED SYSTEM 

The purpose of this section is to derive the equations needed for our numerical calcula- 
tions of the optical conductivity. In linear response theory, the real part of the conductivity 
tensor a^ v is related to the imaginary part of the current-current correlation function 
by§ 

Re[<v(u;)] = -Im[n^(cu)] (1) 

UJ 

where the Greek subscripts refer to spatial components and uj is the frequency. The factor of 
2 is to account for electron spins, which will otherwise be ignored for the rest of this article. 
In the superconducting state, the spins are dealt with explicitly in the calculation of II (a>). 
In a previous article we evaluated II for a two-layered tight binding system0. The intention 
of our earlier calculations was to find the penetration depth, so that II was evaluated at 
zero frequency (uj = 0) and in the superconducting state. Here we will evaluate II(a;) at 
finite frequency and in the normal state. As before, the model is a two-layer tight binding 
model, so that the calculation is very similar to our earlier one. For this reason, the reader 
is referred to our earlier work for details of the calculation which are not shown here. 
In previous work, we showed that@ 

II F (w n ) = e 2 ^ 7T H Tr ( G ( k ' iuJ rn - iVn) 
P m ^ k 

x 7 M (k,k)G(k;^ m ) 7 ,(k,k)] (2) 

where v n and u) m are the boson and fermion Matsubara frequencies respectively, G(k; iu m ) 
are the thermal Green's functions and 7^ are the electromagnetic vertex functions. This 
result is essentially the standard result,@ with the exception that here the Green's functions 



4 



and vertex functions are 2x2 matrices whose diagonal elements (eg. Gn) describe properties 
of a single layer and whose off-diagonal elements describe the interlayer coupling. The trace 
in Eq. (Q) is over the matrix product contained in the square brackets. In order to find the 
optical conductivity, we need to find explicit forms for G and 7^, and we begin by introducing 
our model for the two layer system. 

The model we are going to present describes a metallic system with two types of layer 



stacked along the c-axis (or equivalently, the z-axis). In section |TT| one of the layers is a two- 
dimensional plane layer, while the other is a one- dimensional chain layer. In section [V] both 
layers are plane layers. We define the operators cik and C2k to be the annihilation operators 
for the two types of layer. The wavevectors k are three-dimensional. The HamiltonianllUHll 
for our model is 



c lk c 2k 



h(k) 



Clk 



C 2 k 



(3) 



with 



h(k) 



6 t 
t* 6 



(4) 



and where £i(k) and &{k) are the energy dispersions for the two types of layer, and t(k) 
connects the layers through single electron hopping. The specific forms of £1, £2 arid t will 
be given in sections [TT1| and |V[ 

The Hamiltonian matrix h(k) is diagonalised by the unitary matrices U(k), so that 



e + (k) 
c_(k) 



lP(k)h(k)U(k), 



(5) 



where 



U(k) 



(6) 



and where e±(k) are the eigenvalues of h(k) (ie. e±(k) are the band energies), 
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6 + 6 , 
e ± = — t; — ± 



\ 



6-6 



t 2 . 



(7) 



The Fermi surfaces are the solutions of the equations e±(k) = 0. As an example, one possible 
Fermi surface for the chain-plane system — which is discussed in more detail in section [HI] — is 
shown in Fig. [I]. 

The single particle Green's function is determined from the Hamiltonian: 



G(k; iuj ri 



(8) 



iu n -6 -t 

-t* iu n - £ 2 

If we wish to include impurity scattering, then the simplest approach is to introduce a 
scattering rate Y 



(9) 



iuj n - 6 + «rsgn(u; n ) -t 

-t* iu n - £ 2 + zrsgn(uy 

We will assume that T is independent of frequency and momentum, but that it varies linearly 
with temperature as is observed in the copper oxides for the in-plane conductivity. At T = 0, 
the system is in the clean limit (r = 0) and at T = 100 K, T = 10 meV. The scattering rate 
is related to the quasiparticle lifetime by T = h/2r. 

Finally, we need to find the vertex function 7 At (k, k). In previous workSH it has been 
shown that in the tight binding model, 7 /1 (k, k) is just the gradient of the Hamiltonian 
matrix in k space: 



7,(k,k) = i- 6 ' 



h dk,. 



6 t 
t* 6 



(10) 



Equations (0), (|9]) and (|T0| ) are sufficient to calculate the optical conductivity. In their 
current form, however, they are not very revealing and it is difficult to understand the results 
of our numerical calculations without some further work. We will make two manipulations 
in order to make the formula for the conductivity more transparent. The first is to write 
the conductivity in terms of spectral functions, instead of Green's functions: 
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g 2 ^ „oo 

Re[o-^(w)] = y rfx Tr [A(k; x)7 M (k, k) 



f{x)-f{x + hu) 



(ii) 



where the spectral function A(k;u) is defined by analytically continuing G(k,iu n ) to the 
real axis: 

A(k; w) = i [G(k; + iO) - G(k; uj - iO)] 

and where f(x) is the Fermi function. 

The second manipulation is to make a change of basis, so that the eigenstates of the 
Hamiltonian, rather than the eigenstates of the isolated layers, are used as the basis states. 
In other words, we will evaluate 

e 2 h 



2nQ 



E 



dx Tr 



x A(k; x + hu)%(k, k) 



A(k;x)%(k, k) 
f(x) - f(x + huj) 



(12) 



where we introduce the notation throughout this work 0(k) = (k)0(k)U(k), with 0(k) 
a 2 x 2 matrix. 

Equation (|TJ) is more useful than Eq. (0) for two reasons. The first is that the spectral 
function has a simple interpretation as the density of states. The second is that A has a 
particularly simple form. The price we pay is that the vertex function 7 M is more complicated 
to evaluate than 7 M . 

The vertex function 7 M can be found by explicitly performing the matrix multiplication 
Z7T(k)7^Z7(k), and we will explore it in detail in the following sections. The spectral function 
can be evaluated easily here. In the new basis 

G(k]iu n y 1 = 
iu n - e + + iTsgn(u n ) 

iu n — e_ + «rsgn(o; n ) 



(13) 
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and 



A(k;uj) 



2r 

( w _ e+ )2 +r 2 







2r 



(14) 



(w-e_) 2 +r 2 ' 

The diagonal elements of the spectral function can be interpreted as the density of electronic 
states in the bands. 



III. PLANE-CHAIN MODEL 

In this section we describe a simple model containing a plane layer and a chain layer. 
The model is meant to incorporate the most important features of the chain-plane coupling 
in YBC07_,5. In the discussion which follows, it will be made clear that the interband terms 
in the optical conductivity (which are important for a zz ) are sensitive to the specifics of the 
band structure, which we cannot hope to describe correctly with our model. However, we 
will still able to draw a number of general conclusions which should apply to models with a 
more accurate description of the unit cell. 

For this work we will take the dispersion in the plane layer to be 

£1 = —2ti[cos(k x a) + cos(k y a) 

—IB cos(k x a) cos(k y a)} — /ii, (15a) 

and the dispersion in the chain layer to be 

£ 2 = -2t 2 cos(k y a) - n 2 , (15b) 

where a is the lattice constant in the a (or equivalently x) and b (or y) directions. There 
have been a number of attempts to fit the Fermi surface of the plane layer either to band 
structure calculations^^! or to angle resolved photoemission experiments^, however, we 
are unaware of any work that attempts to find a phenomenological chain Fermi surface. In 
any case, a more accurate Fermi surface will not affect our basic (qualitative) conclusions, 
although the quantitative results of our calculations are quite sensitive to the band structure. 
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The simplest model of interlayer coupling is 

t = -2i ± cos(M/2), (15c) 

where d is the lattice constant in the c (or z) direction. The form of Eq. ( |15cj ) makes the 
assumption that we have coherent single electron transport along the c-axis. 

There have been many models of incoherent transport between layer sJIBBii'i A review 
of some of the models has been provided by Cooper and Gray! As these models are not 
directly relevant to our present work, it will be sufficient here to provide a few highlights. In 
the work of Graf et al.,@ no contribution from coherent transport in envisaged (i.e. t± = 0) 
and the c-axis transport proceeds entirely through incoherent elastic scattering. In this 
case, the resistivity along c- is inversely proportional to that in the afe-plane. On the other 
hand, Kumar and JayannavarB envisage tunnelling between layers but in the limit that 
the tunnelling time is larger than the in-plane scattering time. The effective transverse 
tunnelling matrix element is modulated and reduced by the in-plane scattering and the c- 
axis conductivity becomes proportional to the a — b plane scattering time. Leggetti envisages 
thermal depairing between layers for the case t± < ksT while Anderson and Zhoui consider 
the possibility that the Cu02 planes cannot be described in Fermi liquid theory and involve 
spin and charge separation. This implies blocking of c-axis transport. Finally, Rojo and 
Levin! consider the possibility that c-axis transport proceeds through the combination of t± 
and incoherent transport due to elastic (impurity) and inelastic (phonon) assisted processes. 

Here only the coherent contribution is included. It should also be remembered that we 
have oversimplified the internal structure of the unit cell in YBCO, particularly since we 
have only included one type of hopping mechanism which is paramaterised by t ± . In reality, 
there could be several types of hopping process occuring. For example, the coupling between 
planes and chains could be coherent while the intercell coupling could be incoherent. These 
issues are not treated here and go beyond the scope of our work. They are mentioned, 
however, so that the reader understands clearly the limitations of our work. We do feel, 
however, that near optimal doping YBCO is one of the few materials to display coherent 



transport between layers, and that our description of the interlayer coupling is reasonable. 
The vertex functions are 



7x(k, k) 



7y(k,k) 



Vlx 



Vly 
V 2 y 



and 



7,(k,k) 



v± 
v ± 

where = h~ 1 d£i/dk fM , and v±_ = K~ l dt/ dk z . 
The change of basis is simple to perform and 



7r(k,k) 



e i - e_ 



%(k,k) 



Vly 



+ 



V 2 y 



6 - e_ t 
t e+ - 6 

6 - e- * 
* e+ - 6 

6 - e_ t 



and 



7z(k,k) 



(16a) 



(16b) 



(16c) 



(17a) 



(17b) 



(17c) 



2t 6 - 6 
6 - 6 2t 

In the derivation of Eqs. (|T7|), use was made of the fact that t(k) > 0. We are now able to 
evaluate the real part of the conductivity tensor [Eq. (|T2D1 using Eqs. ( PI) and (|T7j). 

In Fig. the conductivity is plotted as a function of frequency for a number of different 
temperatures. The strength of the interlayer coupling is relatively weak (t± = 5 meV) and 
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this is reflected in the small magnitude of the c-axis conductivity. We point out that the 
case presented in Fig. |2| provides a qualitative description of what is observed underdoped 
YBCO6.7-00 To begin with, the conductivity in the a direction has a traditional Drude- 
like structure. Since the scattering rate scales linearly with temperature, the width of the 
Drude-peak also scales linearly with temperature, and at T = 0, the conductivity is a 5- 
function at uo = 0. The conductivity in the b direction also has a predominantly Drude-like 
structure but its magnitude is roughly twice that of a xx . This is because there are two 
current carrying channels in the b direction (the chains and the planes) and only one in the 
a direction (the planes). At low temperatures there is a small non-Drude contribution due 
to interband transitions which becomes visible, and at T = 0, the interband transitions are 
the only mechanism for conductivity at finite frequency. 

In contrast with a xx {oj) and a yy (u>), the conductivity in the c direction has a decidedly 
non-Drude behaviour and, instead, appears as a broad background. As we shall see, this 
is entirely due to interband transitions between the plane and chain layers. A further, 
interesting feature in <j zz (uj) is that at high temperatures there is very little structure as 
a function of frequency u, but that structure appears at low temperatures. This can be 
attributed to the linear decrease in the scattering rate with temperature assumed in our 
work; at high temperature the large scattering rate smears out the structure in a zz [uS). 
Perhaps the most striking feature in a zz is the development of a gap as the temperature is 
lowered. This gap is just the gap between the two bands e + and e_, and at high temperatures 
it is filled in by the large scattering rate V. 

In Figs. |3| and £|, the optical conductivity is plotted for larger values of the chain-plane 
coupling t±. The larger coupling has almost no effect on a xx . On the other hand, the larger 
coupling enhances the role of interband transitions in a yy , and also increases the Drude (or 
intraband) contribution to a zz . Note that the Drude peak at low lu is not seen in a zz in Fig. 
^| where t± is small, while in Fig. [|, where t± is much larger, the Drude peak is seen as a 
sharp turn upward at low uj. 

The conductivities shown in Figs. [| and £| are not quite appropriate to describe optimally 
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doped YBCO even though the magnitudes of the conductivities are approximately right. 
This is because there is still a gap in the conductivity at low temperatures. However, by 
shifting the chain Fermi surface slightly, as shown in Fig. |5], we are able to eliminate the 
band gap. It is very plausible that a change in doping will change the structure of the chain 
and plane Fermi surfaces so that they cross. In Figs. ^] and the conductivity is shown 
at T = 100 K and at T = K for tj_=10 meV and 20 meV, respectively. In Fig. || the 
conductivity has a Drude-like appearance in both the a and b directions, even though there 
is a large interband contribution to a yy . The interband contribution can be resolved in 
a yy if the temperature (and therefore the scattering rate) is reduced. Unfortunately, this is 
difficult to do experimentally because of the onset of superconductivity near 100 K. In both 
Fig. |^ and |7|, the c-axis response shows a Drude-like peak at low u and high temperature 
(solid curve). On closer examination, the response is a combination of a Drude and flat 
interband contribution. As the temperature is lowered, this second contribution remains 
(dotted curve) and a pseudogap at low lo is resolved (i.e. the conductivity is depressed but 
goes strictly to zero only at u;=0). 



IV. CLEAN LIMIT 

The purpose of this section is to understand the role of interband transitions a little 
better. We will examine why it is that the interband transitions are most important in 
the c-axis conductivity for small interlayer coupling, while they are most important for 
the 6-axis conductivity for large interlayer coupling, and are never important for the a- 
axis conductivity. Finally, we shall derive expressions for the intraband and interband 
conductivities in the clean (r = 0) limit and show how a zz (u) can be interpreted as a probe 
of the band structure. 

Let us start by writing out Eq. (|T2"D, the formula for the conductivity, explicitly: 

Re[(7 w (o;)] = — — / dx A u (k;x)A u (k;x + hu)[y 11 ] 2 + A 22 Ck; x)A 22 (k; x + hu)[% 2 }l 

Z7TS6 , J —oo L 



k 
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+ (in(k; x)i 22 (k; x + hto) + i 22 (k; x)A u {k; x + hu)) [712U721 



f(x) - f{x + fig) 



The first two terms in the integrand describe intraband processes and yield Drude-like be- 
haviour. This fact can be made clear if we remark that 

dU{k) 



[7 ^ = t^: + ^-^ 



f/t (k) . 



(19) 



(where we use the notation ei = e + and e 2 = e_) so that the first two terms in the conduc- 
tivity are 

, , e 2 h 



2lltt 



X 



EE 

± k 

2T 



POO 

/ dx 


~lde± 


1 — 00 


h dk^ 



2T 



x 



(x - e±) 2 + P [hu + x - e±) 2 + r 2 
f(x) - f(x + hw) 



(20) 



In the simple case of a single band system with a spherical Fermi surface, this becomes the 
usual expression for the optical conductivity 



<t(u) = 2e 2 N(0) 



T 



3 1 + cj 2 t 2 ' 



where T = h/2r. 

In the clean limit (r — > 0) 



X" 



+ r 2 



7TS(x). 



and the intraband conductivity becomes 

, , 2ire 2 h 



0"Drude(^) 



EE 

± k 



lde± 



h dk., 



df(e ± ) 
de ± 



(21) 



The Drude peak in the conductivity is now a 5-function centred at u = 0. 

The remaining terms in Eq. ([18]) are the interband terms. In the clean limit the interband 
contribution is (for uo > 0) 



2ne 2 h^ 2 /(e_) - f(e+) 

<wM = — E bi 2 ], — ^ — 

k "r 

xd(huj — e + + e_). 
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(22) 



In order to get some sense of the relative magnitudes of 0"Drudc(^) and center (<^) we 
compare the magnitudes of [h~ 1 de±/dk fl ] 2 and [712]^- In the analysis which follows, we will 
compare these two quantities for a xx , a yy and a zz . 

In the x-direction, 



lde 



± 



and 



Via 



6 - e T 

' e + — e_ 



Via 



e + — e_ 



Throughout most of the Brillouin zone, |£i — £ 2 | 3> 2\t\ and 

4-2 

e + ~ max^i,^) + 



(6 - 6) 



min(£i,£ 2 ) - 



(£i-6) 2 

from which it follows that the largest contribution to ODrudc(^) will be of the order v\ x x 
[1 + 0(t 2 /(^i — C2) 2 )], where the correction term is of the order 1/25 throughout most of 
the Brillouin zone. In the same way, we can immediately see that throughout most of 
the Brillouin zone, the interband contribution to a xx will be v\ x x 0(t 2 /(^i — £ 2 ) 2 )- Of 
course, in the region of the Brillouin zone where the chain and plane Fermi surfaces are 
close together, the above argument does not hold. However, the value of v ix is small in 
this region of the Brillouin zone and we conclude that interband processes do not make a 
significant contribution to cr xx . 

Much of the above argument holds for a yy as well, and we can conclude that throughout 
most of the Brillouin zone, interband processes are unimportant. Unlike the case of a xx , 
however, v± y and v-i y are not small in regions of the Brillouin zone where interband transitions 
are significant, and there will therefore be small but noticeable non-Drude contribution to 
the 6-axis conductivity. 

The situation is somewhat reversed for the c axis conductivity, where we can show that 
interband processes play a dominant role. In the z-direction 
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~lde± 

k dk 7 



2t 



v±- 



and 



6-6 



Inspection of these two equations shows that while the interband contribution still scales 
as t\ (through the factor of v±), the Drude component now scales as t\. In other words, 
Drude contribution to the conductivity is smaller than the interband contribution by a factor 

t 2 /(6-6) 2 . 

In summary, then, the conductivities scale with t± as follows: For a xx and a yy , we find 
that the Drude part scales as (t±)° and the interband part scales as t\, while for a zz the 
interband part scales as t\ and the Drude part as t\. 

We will finish this section with a brief mention of the usefulness of probe of the 

band structure. The factor 5(huj — e + + e_) in Eq. fl22|) means that the interband conductivity 
is a probe of the joint density of states of the two bands. The sum in Eq. ( p2| ) is weighted 
by the thermal factor /(e_) — f{e + ), which restricts the transitions to be between filled 
and empty states. At zero temperature, this term restricts the sum to lie in regions of the 
Brillouin zone where e + e_ < (ie. e_(k) is a filled state and e + (k) is an empty state). In 
Fig. m this means that the interband transitions occur in the area contained by the two 
Fermi surface curves. The constant energy-difference contours are also shown in Fig. [l] as 
the dotted lines. 

When we compare the chain-plane system with the bilayer model in the next section, we 
will see that the chain-plane system is somewhat special in that the joint density of states 
is spread over a broad range of energies. In simple terms, the interband contribution to the 
c-axis conductivity exists over an energy scale of 1 eV because the chain and plane bands 
have such different structure, and the energy difference e + (k) — e_(k) takes on all different 
values on the energy scale of an eV in the Brillouin zone. In the bilayer system, however, 
where the two layers are identical, the energy difference is restricted to a narrow range of 
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values, and the interband contribution to the optical conductivity will result in a narrow 
peak. 



V. BILAYER MODEL 

In this model, we have two planes per unit cell. The planes have equivalent band struc- 
tures, but are made inequivalent by their spacing. In other words, the planes are alternately 
spaced by distances d\ and d 2 along the c-axis, where d\ + d 2 = d is the unit cell length. In 
the special case d\ = d 2 , the model reduces to a single band model as we shall show below. 

The energy dispersion in the planes is £1 = £ 2 = £. For our numerical calculations, we 
will take 

£ = — 2ti[cos(k x a) + cos(k y a) 

—IB cos(k x a) cos^a)] — (23a) 

In this section, we will also calculate the conductivity analytically using the simpler band 
structure 

e = (23b) 



where vf is the Fermi velocity at the (circular) Fermi surface, k\\ = yjk% + ky and kp is the 
value of at the Fermi surface. 

The interlayer coupling term takes the formll! 



t(k) = t ±1 e ik * dl + t ±2 e- lk * d2 , (23c) 

where we expect that if d\ < d 2 , then t±i > t± 2 - If d\ = d 2 and t\ = t 2 , then t(k) reduces to 
Eq. (|15c|). The band energies are 



e ± = C(k)±|t(k)|, (24) 



and the Fermi surface is shown in Fig. H for the case where £ is given by Eq. (|23a|). In the 



bi-planar model, the bands are split by 2|t(k)| so that the maximum band energy difference 
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is 2 1 + tj_2 1 and the minimum band energy difference is 2|tj_i — t± 2 \- Since t±\ and t± 2 
are typically much smaller than the bandwidths, this means that the interband contribution 
to o zz will be over a small range of energies. As we will see from our numerical work, the 
bilayer is not a likely source for the broad, experimentally observed, c-axis response. 
The unitary matrix which diagonalises the Hamiltonian is now 



U(k) 



V2\t\ 



-t -t 
-\t\ \t\ 



(25) 



The vertex function for the in-plane conductivity a xx is 





7 x (k,k) = 7a,(k,k) 



v x 



(26) 



where v x = h d£/dk x . We can see that there will be no interband contribution to a xx (u) 
since j x has no off-diagonal matrix elements. In fact, Eq. ( |T2"D becomes 

e 2 h 



Re[a xx (u)] 



2ttO 



EE 

± k 

2r 



- 2 f(x) - f(x + huj) 

ax v 

-oo nu) 



2T 



x- 



x - e±) 2 + r 2 (x + huj - e ± ) 2 + r 2 ' 



If we take £ from Eq. 



then we get, in the usual way, the Drude conductivity 



Re[a x 



4e 2 ftiV|| v 2 F 



T 



(27) 



(28) 



d 2 u 2 t 2 + r 

where N\\ is the two dimensional density of states for a single layer. In three dimensions, with 
a cylindrical Fermi surface, the density of states for the bilayer is iV = 2N»/d. Numerical 
calculations of the optical conductivity, shown in Fig. [9|, are in qualitative agreement with 
Eq. ©. 

The vertex function for the conductivity along the c-axis is 

v ± 
v*, 



7 2 (k,k) = ^ (k) 



U(k) 



Re(v z t*) -i lm(v z t*) 
i Im(v z t*) - Re(v z t*) 



(29) 
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where v z = h 1 dt/dk z . The diagonal elements in j z are 

Re(v z t*) 1 d\t\ 



t 



Ti dk z 

1 t±ltj_2 

h \t\ 



d sin(k z d) 



< - min(t_Li, t ±2 ) d sin(fc z d), 

so that the intraband (or Drude) conductivity is limited by the lessor of t±i and t± 2 . In 
other words the Drude current is limited by the weak link along the c-axis. The off-diagonal 
elements in are 

lm(v z t*) 1 



Ti 



[t ±1 di — t ± 2^2 



+t±it± 2 (di - d 2 ) cos(k z d)]/\t\ 



(30) 



which vanishes when t\ = t 2 and d\ = d 2 - Unlike the intraband conductivity, the interband 
conductivity does not become small as either t\_\ or t\_ 2 vanishes. In other words, the 
interband contribution to o zz persists even in the limit that the bilayers become isolated 
from their neighbours. 

The conductivity along the c-axis is 



Re[a zz (u)] 



e 2 h 
2nQ 



E 



tlx f{x) ~ f{x ±M(^ 



2T 



2T 



[x - e±) 2 + T 2 (x + hcu- e±) 2 + T 2 



Re\v z f 



t 



+ E 



2T 



2T 



lm\v,t* 



1 2~ 



± (x- e±) 2 + T 2 (x + huj - e T ) 2 + T 2 [ \t\ 
The first term in the curly braces gives the intraband conductivity while the second gives 
the interband conductivity. We can proceed further if we take £ to be of the form given in 
Eq. ( |23b| ). Then we find that 



(31) 



Re[a zz (u)} = < 2 



T 



(jJ 2 T 2 + 1 



E 



T 



\t\~ 



, ± x T 2 (LU±2\t\/h) 2 + l 

where ()k z denotes an average over k z . 

The first term in Eq. ( |3^ ) has the usual Drude frequency dependence, weighted by an 
average Fermi velocity. The second term is the interband term. Its frequency dependence is 



(32) 



an average over k z of Lorentzians centred at hu = 2\t(k z )\. This is what is seen in Fig. (^j). 
In the limit that we have a single bilayer (t± 2 - > 0), |t(& z )| — > t_u and 

^ r , 2e 2 N ll tl l d 2 l t . . 

Re | Mw)] _ _____ (33) 

The c-axis response becomes a Lorentzian centred at hu = 2t±i. This case has been studied 
in detail by Gartstein, Rice and van der Marel.0 



VI. SUM RULES 

In this section we will discuss the partial sum rule for the conductivity within our model. 
The full sum rule 



irne 



2m 



duo <j(uj ) 



(34) 



where n is the electron density and m is the bare mass, is well known. Equation Q53|) is 
often used as a means of measuring the electron density in the cuprate superconductors. 

In any practical evaluation of Eq. (|34|), it is necessary to impose a cutoff in the frequency 
integration, and often this cutoff is taken to be below the onset of interband transitions 
in the a and b directions, so that only the Drude-like response is counted. If we define, 
therefore, a plasma frequency tensor by 

,2 



^ = / duj cr Drude (u) 
o JO 



and consider r — > for simplicity then, from Eq. (f2l|), 

1 de± de± 



(35) 



ire 



8 n 

Integrating by parts in gives 



EE 

± k L 



" ± k 



h 2 dk^ dk p 



1 d 2 e ± 
h 2 dk a dK 



df(e±) 
de ± 



(36) 



/(6±)- 



(37) 



Now, we can define an average effective mass tensor for each of the bands by 
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where 



so that 



n± il X 



1 d 2 e± 



k 2 dkpdK 



(38) 



k 



(39) 



CO 



■ptiv 



(40) 



We emphasize here that M±\ v is an average of the effective mass tensor over all filled states, 



and is not just the effective mass tensor at the Fermi surface. As a result of this, M±p„ 
depends on the filling of the bands. For this reason, we suggest that it may actually be 
difficult to determine changes in the electron density with doping in the high T c materials. 



VII. CONCLUSION AND DISCUSSION 

We have calculated the optical conductivity for a simple layered model of YBC0 7 _5, in 
which each unit cell consists of a two-dimensional plane layer and a one-dimensional chain 
layer. The model contains two important pieces of physics. First, we assumed that the layers 
are coherently coupled and measured the strength of the coupling with the parameter t±. 
As we mentioned earlier, this is different from the more usual point of view which ascribes 
the broad background and absence of a Drude peak in the c-axis conductivity to incoherent 
c-axis transport. The second important piece of physics is the scattering rate T, which was 
assumed to vary linearly with T. 

In a multilayer system, we found that there are two contributions to the optical conduc- 
tivity: intraband and interband. For the in-plane conductivity (a xx and a yy ) we found that 
the intraband conductivity dominates the response, resulting in Drude-like conductivities. 
In fact, we showed that the intraband terms scale as (t±)°, while the interband terms scale 
as t\ for small t±. On the other hand, we found that the c-axis response is dominated by the 
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interband contribution (which still scales as t 2 ± ) since the intraband contribution scales as t\. 
For this reason, we found that the c-axis optical conductivity has a non-Drude appearance. 

In the plane-chain model, the a zz consists of a broad featureless response at high tem- 
peratures. As the temperature is lowered, structure develops in the conductivity. In some 
cases, we found that a pseudogap-like structure appeared. The broad range of frequencies 
spanned by a zz is the result of the plane and chain bands having different geometries so that 
transitions between the two bands span a wide range of energies. In contrast, we showed 
that a bilayer system in which the two layers have identical band structures results in a 
c-axis conductivity which is finite only over a narrow range of frequencies. 

The other feature of a zz — that structure appears as T is lowered — is the result of having 
a temperature dependent scattering rate. At T = 100 K, for example, structure in a zz is 
smeared out over 20 meV, while at T = 10 K, structure can be resolved on a scale of 2 meV. 

Finally, with this model, we have been able to comment on changes in a zz with doping 
over the range YBC0 6 . 7 to YBC0 7 . We have suggested that in slightly underdoped YBCO, 
there is a band gap between the plane and chain layers, which is reflected in a gap in the 
interband conductivity. Above a certain temperature, however, the band gap is hidden by 
the large scattering rate which smears out quasiparticle energies by more than the band 
gap. As the temperature is reduced, the quasiparticle energies become better defined and 
the band gap appears in the conductivity. The gap is similar in its appearance to the 
pseudogap observed in YBC0 6 . 7 . At optimal doping, we have suggested that two changes 
must be made to the model. The first is that the strength of the chain-plane coupling must 
be increased — thus increasing the importance of the Drude contribution to a zz . The second 
is that the chain and plane Fermi surfaces must be made to cross, eliminating the pseudogap. 

We point out that the exact shape of the interband contribution as a function of energy 
was found to be quite sensitive to details of the Fermi surfaces involved. This means that 
c-axis conductivity measurements could, in principle, be used to get information on the 
energy bands as well as on the filling factors for chains and planes and on their changes with 
oxygen doping. 
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While the sensitivity of the frequency dependence of cj c {uj) to band structure could, 
in principle, be used to get spectroscopic information on electronic structure, it should 
be emphasized that the band structure used in our work is grossly oversimplified and so 
some of our detailed predictions cannot be applied directly to YBa 2 Cu 3 O x . The qualitative 
features obtained and emphasized in this conclusion are, however, expected to be robust and 
remain in more complex models. Such calculations will need to employ more realistic band 
structures and, perhaps more importantly, consider the issue of intercell coupling which 
could be incoherent and quite different from the chain-plane hopping. Nevertheless, our 
model does exhibit many of the features observed in experiments on YBCO at optimum 
doping as well as underdoped and overdoped cases. 
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FIGURES 



FIG. 1. Fermi surface for the chain-plane model. The Fermi surface is shown for k z = (solid 
curve) and k z = ir/d (dashed curve). When k z = ir/d, the chain-plane coupling i(k) vanishes and 
the Fermi surface is that of the isolated chain and plane layers. The dotted curves are lines of the 
constant energy difference u = e+ — e_. The band structure parameters chosen for this case are 
{h,t 2 , m, /J,2, t±} = {70, 100, -65, -175, 20} meV and B = 0.45. 

FIG. 2. Normal state optical conductivity in the (a) x (b) y and (c) z directions for the 
plane-chain model. The conductivity is shown for temperatures 200 K (solid curve), 100 K 
(dot-dashed curve), 10 K (dashed curve) and K (dotted curve). The scattering rate is 1/r = 20 
meV at 100 K and it scales linearly with temperature. The conductivity in the z-direction is 
dominated by interband processes and has a non-Drude appearance, while the conductivity in the 
x and y directions has a predominantly Drude-like behaviour, although there is a small interband 
contribution to a yy . The band structure parameters are the same as in Fig. [|, except that t± = 5 
meV. 

FIG. 3. The normal state conductivity is shown for the (a) x, (b) y and (c) z directions at 
100 K (solid curve) and K (dotted curve). The band structure differs from that of Fig. ^ by the 
magnitude of the chain-plane coupling, which is t± = 10 meV here. 

FIG. 4. The normal state conductivity is shown, as in Fig. |, but with t± = 20 meV. The 
conductivity in the z direction has a significant Drude part, while the conductivity in the y direction 
has a significant interband part. 
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FIG. 5. The Fermi surface is shown for the chain-plane model for a case where the chain 
and plane Fermi surfaces (dashed curves) cross. This case is qualitatively different from the 
case shown in Fig. |l| because there is no band gap. As a result, a pseudogap is not expected 
in the c-axis optical conductivity. The band structure parameters chosen for this case are 
{h,t 2 ,m,lJ,2,t±} = {70,100,-65,-130,20} meV and B = 0.45. 

FIG. 6. The optical conductivity is shown along the (a) x, (b) y and (c) z directions at T = 100 
K and at T = K. The T = K optical conductivity is entirely due to interband transitions at 
finite frequencies. In a xx there is only a very small interband contribution, while in a yy and a zz , 
the interband contributions are substantial. At T = 100 K, however, it is difficult to distinguish 
the interband contribution from the Drude contribution because of the large scattering rate. There 
is also a large Drude contribution to a zz and, and there is no pseudogap, although the interband 
conductivity still falls to zero linearly with lo. The band parameters are the same as in Fig. ||, 
except that t± = 10 meV. 

FIG. 7. The optical conductivity is shown, as in Fig. H but with t± = 20 meV. 

FIG. 8. The Fermi surface is shown for the bilayer system. In this model, there are two plane 
layers per unit cell. The dispersions in the layers are the same, but they are spaced distances d\ and 
c?2 apart in alternating fashion. The Fermi surfaces are shown for k z = (solid line) and k z = ir/d 
(dashed line). The minimum and maximum coupling strengths are 2\t±i — t±2\ and 2|tj_i + t±2\ 
respectively. The model parameters are {t\, /ii, t±i, t±%} = {70,-65,20,10} meV, d\ = 0.3d and 
B = 0.45. 
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FIG. 9. The optical conductivity is shown for the bilayer system (a) in the a and 6-directions and 
(b) in the c-direction. The conductivity is shown at T = 100 K (solid curve) and at T = K (dotted 
curve). The in-plane conductivity has a Drude shape, while the c-axis conductivity has a large 
interband contribution. At T = K, the in-plane conductivity vanishes at finite frequency. Along 
the c-axis, however, the interband contribution remains for frequencies 2|tj_i — tj_2| < Tiu)2\t±\+t±_2\- 
The model parameters are the same as those in Fig. @. 
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